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Abstract —This paper proposes a belief propagation (BP) mes¬ 
sage passing algorithm based joint multiple symbol differential 
detection (MSDD) and channel decoding scheme for noncoherent 
differential ultra-wideband impulse radio (UWB-IR) systems. 
MSDD is an effective means to improve the performance of non¬ 
coherent differential UWB-IR systems. To optimize the overall 
detection and decoding performance, in this paper, we propose 
a novel soft-in soft-out (SISO) MSDD scheme and its integration 
with SISO channel decoding for noncoherent differential UWB- 
IR. we first propose a new auto-correlation receiver (AcR) 
architecture to sample the received UWB-IR signal. The proposed 
AcR can exploit the dependencies (imposed by the differential 
modulation) among data symbols throughout the whole packet. 
The signal probabilistic model has a hidden Markov chain 
structure. We use a factor graph to represent this hidden Markov 
chain. Then, we apply BP message passing algorithm on the 
factor graph to develop a SISO MSDD scheme, which has better 
performance than the previous MSDD scheme and is easy to be 
integrated with SISO channel decoding to form a joint MSDD 
and channel decoding scheme. Simulation results indicate the 
performance advantages of our MSDD scheme and joint MSDD 
and channel decoding scheme. 

Index Terms —Ulta-wideband (UWB), multiple symbol differ¬ 
ential detection (MSDD), channel decoding, belief propagation 
(BP), factor graph. 


I. Introduction 

U LTRA-wideband impulse radio (UWB-IR) is served as a 
promising candidate for location-aware indoor commu¬ 
nications, wireless sensor networks and wireless personal area 
networks. Previously, UWB-IR earned significant attentions 
in both academia and industry |Tj. However, the implementa¬ 
tion of optimal coherent receiver for UWB-IR systems faces 
many challenges. UWB channels usually contain hundreds 
of multipath, due to the rich scattering indoor environments. 
The optimal coherent receiver required to capture multipath 
energy is the famous Rake receiver Since the UWB 
channel is characterized by the dense multipath, we need a 
large number of Rake fingers to capture a significant part 
of the signal energy 0. The implementation of so many 
Rake fingers and the associated channel estimation on the 
corresponding multipaths involve intensive complexities 0. 
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Moreover, such Rake receiver is very sensitive to timing-jitter 
|5). These challenges make it difficult and costly to realize the 
optimal coherent receiver for UWB-IR systems. 

To obviate the complicated treatments on UWB channels, 
the suboptimal noncoherent receivers are proposed (6). The 
typical noncoherent UWB-IR schemes are differential |7j and 
transmitted-reference (§ UWB-IR systems, both deployed 
with the analog autocorrelation receiver (AcR) that does not 
require Rake receiver and explicit channel estimation. Due 
to their good performance-complexity tradeoff, noncoherent 
receivers are now more popularly used in UWB-IR systems. 
However, they suffer from some performance degradations 
compared with coherent receivers. 

Multiple symbol differential detection (MSDD) is an ef¬ 
fective means to improve the performance of noncoherent 
differential UWB-IR systems. The theoretical framework of 
MSDD is the maximum-likelihood (ML) sequence detection, 
which is firstly introduced to detect a block of differential 
MPSK symbols over additive white Gaussian noise (AWGN) 
channel in j9j- Applying MSDD to differential UWB-IR 
systems is considered in ng-GT Works & G3 consider 
the application of sphere decoding algorithm to fulfill a low 
complexity MSDD for differential UWB-IR systems. With 
the same purpose, work proposed a Viterbi algorithm 
based MSDD scheme and works jTg-og proposed decision- 
feedback MSDD schemes. 

Wireless communication systems are susceptible to various 
impairments, such as noises, interferences and channel fading. 
Channel codes are usually employed to protect the transmitted 
symbols over possible errors. The decoding of most powerful 
channel codes that can approach the Shanon capacity depends 
on iterative algorithm, where iterations are performed between 
soft-in soft-out (SISO) modules [|T7]| . These MSDD schemes 
mentioned earlier, however, are all about to detect the hard 
decisions of the differential modulated UWB-IR signals, which 
is not compatible with SISO channel decoding. Recently, the 
work (18) investigates SISO MSDD for UWB-IR systems and 
incorporate it with SISO channel decoding. 

In this paper, we propose a new SISO MSDD scheme for 
noncoherent differential UWB-IR systems. Even without con¬ 
sidering channel encoding, there are memories introduced by 
the differential modulation to all modulated symbols through¬ 
out the whole packet. In |18j, the SISO MSDD processes sig¬ 
nal samples block-by-block and it just ignores the information 
dependencies among different blocks. This leads to informa¬ 
tion loss. In this paper, in contrast, the proposed SISO MSDD 
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scheme calculates the soft information of one symbol by ex¬ 
ploiting the signal samples from the whole packet. We propose 
an AcR architecture to enable this scheme. The proposed AcR 
architecture correlates the received UWB-IR signals, and does 
not need explicit channel estimation. Thus, the proposed AcR 
is a noncoherent receiver. Moreover, it can exploit the signal 
dependencies among the whole packet. The joint probability 
function of the data symbols and the correlation samples in 
the packet has a hidden Markov chain structure. We use a 
factor graph to represent this hidden Markov chain of the 
signal probabilistic model. We then develop our SISO MSDD 
scheme using the belief propagation (BP) message passing 
algorithm which implements sum-product rule on the factor 
graph G3- The proposed MSDD scheme has better detection 
performance than the previous block-by-block MSDD scheme 
in G3-G3- We also consider the channel decoding for 
noncoherent differential UWB-IR systems. Since BP message 
passing algorithm is also employed as the decoding algorithm 
for many channel codes, we integrate the proposed MSDD and 
the BP messag passing algorithm for SISO channel decoding. 
The outputs of the MSDD are fed to the inputs of the channel 
decoding, and vice versa, in an iterative manner. The main 
contributions of this paper are summarized as the follows. 

1) New noncoherent AcR architecture for generating cor¬ 
relation samples. We propose a new noncoherent AcR 
architecture to sample the received UWB-IR signal. 
Compared with the existing AcR in (10)-|(Tg, the 
proposed AcR can exploit more signal dependencies 
imposed by the differential modulation. The proposed 
AcR results in a hidden Markov chain model for the 
signal of the whole packet. 

2) BP message passing algorithm for SISO MSDD scheme 
(joint MSDD and channel decoding scheme). We ap¬ 
ply the BP message passing algorithm to the fac¬ 
tor graph that represents the hidden-Markov-chain-type 
signal model for deriving SISO MSDD scheme. The 
proposed MSDD scheme is a bidirectional algorithm 
that consists of a forward and a backward message 
passing. Furthermore, we integrate the proposed MSDD 
with the BP messag passing algorithm for SISO channel 
decoding, and we achieve an iterative algorithm for the 
joint MSDD and channel decoding scheme. We believe 
this is the first time that applies BP message passing 
algorithm to noncoherent differential UWB-IR systems. 

3) Performance evaluations by simulations. Simulations are 
performed to validate and evaluate the proposed scheme. 
The performances of uncoded and coded system under 
the environments of UWB multipath channel are evalu¬ 
ated. The results indicate the performance advantage of 
the proposed scheme over other existing schemes. 

The rest of this paper is organized as follows. The system 
model of differential UWB-IR system is described in Section 
II. Section III introduces the joint MSDD and channel de¬ 
coding scheme. Section IV shows simulation results. Finally, 
conclusions are drawn in Section V. 



Fig. 1. A block schematic diagram of the system model. 

II. System Model 

In this section, we present the system model for UWB-IR 
communications. A block schematic diagram of the system 
model is shown in Fig. [I] Adopting binary antipodal pulse 
amplitude modulation (BPAM), the transmitted signal wave¬ 
form is given by 

N 

S (t) = y diUJ s (t - iT s ), (1) 

i =0 

where di £ {±1} is the i th channel symbol, w s ( t ) is the 
symbol waveform with duration T s and N is the packet size. 
We denote the original information bits by bj £ {0,1}, 
j = 1, ■ ■ ■ , K. Through channel encoding, interleaving and 
modulation, these information bits are mapped to the coded 
data symbols a,i £ {±1}, i = 1 The coding rate is 

R = K/N. Finally, the channel symbol di is obtained by 
differential modulation: di = di-icti, i = 1 where 

do = 1 is the reference symbol. UWB-IR transmissions 
usually employ Nf frames to transmit one channel symbol, 
and each frame includes one very short pulse. According 
to this unique aspect of UWB-IR, the unmodulated symbol 
waveform used in (jTJ is expressed as 

Nf - 1 

W s (t) = y: U}(t- jTf - CjT c ), (2) 

3=0 

where co (t) is the ultra-short pulse with the duration T u 
(referred to as the monocycle in literatures), Tf is the frame 
duration and we have T s = NfTf. The sequence { Cj } in 
Q is a user specific time-hopping (TH) code used for the 
purpose of multiple access. Its elements are integers in the 
range 0 < Cj < N c — 1, satisfying Tf > N C T C . T c is the 
duration of an addressable time chip. Since oj (t) has a very 
short duration, T u is typically on the order of nanoseconds, 
the transmitted signal occupies a huge bandwidth. The frame 
duration Tf is usually hundred or thousand times longer than 
If, resulting in a low duty transmission. 

We consider dense multi-path environments, such as the 
industrial and indoor office (20). The channel impulse response 
(CIR) between the transmitter and the receiver is modeled as 
h (t) = aid ( t — Ti ), where S is the Dirac delta function, 

L is the number of resolvable multipath components (MPCs), 
ai and t; is the gain and the delay of the I th MPC, respectively. 

We define the received pulse waveform as: g (t) = oj ( t ) 0) 
h ( t ), where 0 denotes the convolution operator. Then, the 
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received noisy signal waveform is given by 

r (t) = s ( t) 0 h(t) +n ( t ) 

N Nf — 1 

= E d i E g(t-iT a - jTf - CjT c ) +n(t), 

2—0 J =0 

where n (t) is the additive white Gaussian noise process 
with zero mean and a flat two-sided power spectral density 
A ; Vi/2. With T g denoting the maximum delay spread of the 
received pulse waveform g ( t ), the inter-frame interference 
(IFI) is avoided in the received signal (|3]» by letting Tf > 
Tg ■ (.v r i)T c . 


III. Joint Noncoherent Detection and Channel 
Decoding 

A. New Noncoherent Autocorrelation Receiver 

The optimal coherent detection of UWB-IR signals requires 
an implementation of the filter matched to the received pulse 
waveform g (t). However, the complexities of the implementa¬ 
tion of the match filtering and the explicit channel estimation 
constitute obstacles for the practical use of coherent detection 
in UWB-IR systems. Therefore, we focus on noncoherent 
detection that dose not involve the explicit channel estimation 
and the implementation of the match filtering. To improve the 
performance of noncoherent detection, we apply the MSDD 
scheme to the system. UWB channels are quasi-static in 
typical indoor environments [[20]]. This means the channel 
remains invariant over several symbol durations. Relying on 
this feature, MSDD jointly detects a block of M symbols from 
the received signal in the observation window of size M + 1 
symbol durations 0, 0, 01- 

In this section, based on the concept of MSDD, we develop 
an AcR architecture for noncoherent detection of UWB-IR 
signals. We modify the sampling mechanism of the MSDD 
in {TO), 0, 0. Essentially, we still employ the correlation 
principle derived from GLRT criteria in |I2[ to sample the 
received signal; however, we change its sliding mode of 
the observation window. In 0, each time, the observation 
window of size M + 1 will be slid down M symbol durations 
after the current M symbols have been detected. In a different 
mode, we slide the M + 1 size observation window down 
one symbol duration each time. In other words, the current 
observation window overlaps M symbols with the next obser¬ 
vation window. The sampling mechanism of the proposed AcR 
is illustrated in Fig. [2] In the following, we mathematically 
formulate the proposed sampling mechanism, and then explain 
its implications to noncoherent UWB-IR systems. 

From the received signal r (t ) in the i th observation window 
(i — M) T s < t < (i + 1) T s , we obtain the i th sample vector 
Y, = [Y^i, Yi t 2 , ■ ■ ■, Yi t M] T , where Yi, m is the correlation 
sample between the i th and the [i — m) th symbols 


Yi, m = f y(t + iT s ) y(t+(i-m) T s ) dt, 

o (4) 

m = 1,2, ■ ■ ■, M 


The i-th observation window 


Y, _ _.. 


Symbol i-M 


Symbol i-1 

Symbol i 

Symbol i+1 



The i+l-th observation window 

Fig. 2. The illustration for the sampling mechanism of the proposed AcR. 


with the de-spreading signal 


Nf-l 

y(t)= Y r ( t + J T f + CjT c ). 
i =o 


(5) 


After we finish the computation of Y, , the observation window 
is slid down one symbol duration to (i + 1 — M) T s < t < 
(i + 2) T s , from where we will compute the next sample vector 
Yj+i. Since there is no transmission occurring (r (t) = 0 for 
t < 0), we pad some zeros at the rears of the first M— 1 sample 
vectors: Yi tTn = 0 for i = 1, ■ ■ ■, M— 1 and m = i + 1, • • •, M. 
We make some remarks about the proposed AcR to bring out 
its implications. 

1) The AcR receiver can be realized using analog compo¬ 
nents to avoid the ultra high sampling rate in UWB- 
IR systems. The sampling mechanism of the proposed 
AcR produces the samples in Y, by correlating the i th 
symbol with its previous M symbols as expressed in ([Tji; 
then it slides the observation window down one symbol 
duration to produce samples in Y, f ) ; the correlating 
operations are carried out from the first symbol to the 
last symbol of the packet. By this manner, the proposed 
AcR can exploit the dependencies (imposed by the 
differential modulation) among symbols throughout the 
whole packet. 

2) Based on the samples of the whole packet Y,; for 
i = 1, ■ ■ •, N produced by the proposed AcR, we will 
derive the BP message passing algorithm on a factor 
graph for SISO MSDD in the next subsection. We 
employ a factor graph to represent the probabilistic 
model of the samples and the data symbols in the 
whole packet. Then, we apply the BP message passing 
algorithm to the factor graph for SISO MSDD. We also 
combine the BP message passing algorithm for SISO 
channel decoding with the proposed BP message passing 
algorithm for SISO MSDD, resulting in a BP message 
passing algorithm to joint MSDD and channel decoding 
for UWB-IR systems. Iterative message exchange will 
be performed between SISO MSDD and SISO channel 
decoding. A block schematic diagram about this receiver 
structure is also shown in Fig. |T] 

3) Different from ours, the sampling mechanism of 0- 
0, 0 is on a block-by-block basis. The correlation 
operations only try to exploit the dependencies (imposed 
by the differential modulation) among symbols within a 
block of M symbols. However, the symbol dependencies 
between different blocks are ignored. This is a kind of 
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information loss. We will discuss the MSDD scheme 
of |[K)|-[T2|. [181 in Section III.D in detail. Depending 
upon the proposed sampling mechanism in (|2ji where 
blocks overlap some others, the detection of one symbol 
is capable of making use of the information of all the 
symbols throughout the whole packet. Since more in¬ 
formation are collected, it is expected that the proposed 
scheme could have better performance. 


B. The BP Message Passing Algorithm for SISO MSDD 

In this subsection, we derive the BP message passing 
algorithm for performing SISO MSDD using the samples 
delivered from the proposed AcR. Here, we do not consider 
the effect of the channel encoding and we will discuss it in 
Section IH.C. 

Substituting 0-0 into 0 and using the result of differ¬ 
ential demodulation = IIz=i-m+i az ’ we can ex P ress 

the sample Y i rn as 


Yl a z ) N f E g + n i,r 


( 6 ) 


\z=i-m -\-1 


where E g = f* 9 g 2 (t) dt is the captured energy of the 
received pulse and 

T a 


= J y(t + iT s ) n(t + {i- m) T s ) dt 


o 


+ J n(t + iT s ) y (t + (i — m) T s ) dt 
o 
T g 

+ J n (t + iT s ) n(t + (i — to) T s ) dt 
o 


(7) 


is the discrete noise component. It has been shown in [211 


that n^rr, for all i and m can be approximated to mutually 
independent Gaussian random variables with mean zero and 
variance cr^ = NfN 0 E g + WTgNfi/2, where W is the 
bandwidth of the bandpass filter employed at the receiver. 
This approximation is rather well when Nf is large due to the 
central limit theorem. We will also investigate this Gaussian 
approximation using simulations in Section IV. 

It can be concluded from 0 that the signal samples depends 
on the data symbols and the captured energy E g . To obtain the 
knowledge of parameter E g , our receiver employs an energy 
estimation method: 

N M 

yi yi (8) 


i 

Tg — _ 

9 ~zn 2 


j i—t m =1 


where Z is the number of non-zero elements in {Y,}. After 
that, we detect data symbols by substituting E g into the signal 
model ([6]). In the simulation results of Section IV, we will see 
that this straightforward estimation of E g achieves a rather 
good performance. 

Let Y = {Y.;} be the set containing all the sample vectors 
and a = [ai, 02 , • • •, dN] T be the vector of all the data 
symbols. The target of the SISO MSDD is to calculate the 
a posteriori probability (APP) of data symbol a, given Y: 


p(a,i\ Y) oc ^2 p(a)p(Y|a) (9) 

a :~a,j 


for all i, where notation ~o means the summation over all 
data symbols in a except at. The straightforward calculation 
of (|9j> will involve complexity O (2 ;V ), which disastrously 
increases with N. 

For an efficient calculation of 0 we will derive a fac¬ 
tor graph representation for the probabilistic model of the 
system and employ the BP message passing algorithm onto 
the factor graph. To use the BP message passing algorithm 
and factor graphs, we first factorize the globe probability 
function p(Y|a)p(a) into many small local functions. As 
specified by 0, after the signal sampling of the proposed 
AcR, we meet an equivalent discrete memory channel, where 
the i th sample vector Y, depends on the M data symbols 
{ai, cti- 1 , • ■ •, ai-M+i}- Therefore, we can factorize p (Y |a) 
as 


N 

P (Y |a) = [|p(Yi |di, ■ ■ •, di-M+i), (10) 

i= 1 


where 


P i &i— 1) j O'i— M+l ) 


M 

OC Yl ex P 

m= 1 



(riz=i-m+l fl2 ) NfEg 



( 11 ) 

is obtained by using the Gaussian approximation on the 
discrete noise components. Then, we perform a factorization 
on p (a): 


N 

p( a) = n?(«i) 

i=l 

N 

— 0 P di— 1 , * * * 5 di—M +1 1^2—1; di— 2 , 5 di—M ) 

2—1 

( 12 ) 

which intuitively brings out the Markov property introduced 
by the sampling mechanism of the proposed AcR. The first 
equality of ( fl2| ) is due to the independence of di without 
considering the channel encoding] the second equality of (12 1 
is due to 


p(°i) n pm n p(oj) 

p M = —*=1-= - 

n n ?(“j) 

j=i — M j = i — M 

_ P(a-i,a-i-l,"‘ i a i-M + h a i-l, a i-‘2,"‘ , a i-M) 

p(a>i — i,ai— 2 ,-" 

— P (^47 1? ? &i— M+l \Q'i— 1? &i— 2j ? &1—M ) • 

Based on the probability functions © and ( fill ), we 
find that our system is well represented by a hidden 
Markov chain. We now define such hidden Markov chain 
as follows. For the Markov chain, the i th state is S, = 
[as*, di-i, ■ ■ ■, di-M+i] T , the i th input is di and the i th 
output is x, = [di,di-i, - ■ ■ ,di-M+i] T ■ The behavior of 
the Markov chain is defined by the local check functions 
Tj (di, Xj, Sj_!, Sj) that constrains the transitions from state 
* — 1 to state i to be valid. We have I) (a,, Xj, Sj_i, S^) = 
1 when the combinations of its arguments are possible; 
Ti (di, Xj, Sj_i, Sj) = 0 otherwise. The behavior of a Markov 


1 If considering the channel encoding, the independence of a t can be 
approximately achieved by the interleaving operation. 
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chain usually can be illustrated using a trellis chart. In Fig. 
[3] we give an example of such trellis chart that illustrates the 
Markov chain of our system with M = 2. 

Our signal probabilistic model is a hidden Markov chain due 
to that we cannot directly observe the output Xj, we can only 
observe Y, which is a function of x, plus a noise. Hidden 
Markov chains can be represented using factor graphs G3- 
The factor graph representation for the hidden Markov chain 
of our signal probabilistic model is shown in the upper part of 
Fig. 0 where circles are variable nodes for inputs and outputs, 
double circles are variable nodes for state and squares are the 
factor nodes for local check functions. Then, we can apply BP 
message passing algorithm to the factor graph for performing 
SISO MSDD. Since the structure of the factor graph for SISO 
MSDD is a tree, SISO MSDD is exactly implemented by a 
forward message passing and a backward message passing. 
We give the detailed descriptions about these message passing 
operations in the following. 

The above definition for the hidden Markov chain is con¬ 
sistent with the factorization in ( [T2 | i: 

P ) Fj (u;, Xj , Sj_l, Sj) 

P {O'i 1 O'i— 1 j * * * 3 Uj_ M-\-l | &i —i, <Xj —2 , * * * 7 Cti—M ) 


Therefore, we have 


N 


p (a) oc n p(a i )T i (a i ,x l , Sj,Sj_i). (13) 

i =1 

Substituting (JT0|, (13 i into Q leads to 

P(ai\ Y) 


oc 


W_P{ a j) T o ("/• x • B./ i.SjpiY., X, 

3=1 


=«(Si-D 


N 


55 II P( a i) T j ( a H x H S f-i> S f)P( Y fl x f) 

l a i+ i:iv j=i +1 


=P(s 0 

xp(a,i)Ti (a,i, Xj, Sj_i, Si) p (Y, |xj). 

(14) 

As indicated in (14 1 , we can now calculate APPs p(a,i\ Y) 
using a BP message passing algorithm that implements the 
sum-product rule (19) . The BP message passing algorithm for 
calculating p (a,;| Y) is a bidirectional algorithm consisting of 
a forward and a backward message passing. This is similar to 
the BCIR algorithm for decoding convolutional channel codes 


Using the sum-product rule, the forward messsage passing 
aims to recursively calculate messages a (S,) (already defined 
in {Ft}): 

a(Sj) = 

J 2 a (Si-i) P {a-i) Ti (oj, Xj, Sj_i, Sj)p (Yj| Xj) (I 5 ) 


for each i = 1,2, •••,7V, where notaiton )T)_ S denotes the 
summation over all arguments involved in except Sj, p (ai) 
is the a priori information of a*, T) (a,, Xj, Sj_i, Sj) is the 




Fig. 3. The trellis chart for the defined Makov chain with M = 2. The 
digits one each transition represent the associated inputfoutputs): a; (x;). 


previously defined local check function for the state transition, 
and p ( Y,; | x, : ) is the evidence information from observation 
Yj (shown in (fTT}). Similarly, the backward messsage passing 
aims to recursively calculate messages /3 (S,) (already defined 
in {14}): 

f8(Si) = 

^2 2 p Ti^-i , x^r, Sj, Sj-f-i) 

~s l+1 

xp (Yj+i| Xj+i) 

(16) 

for each i = 1.2. ■ ■ ■. N. After the forward and backward 
message passing once in each direction, we can calculate the 
messages that runs fron check node T t to variable node a, 
(denoted by the red solid arrows in Fig. [4} as 

7(ai)= 17) 

22^a (Si-i) /? (Sj) Ti (ai, Xj, Sj_ x , Sj)p (Yj| Xj) (U) 

for each i = 1, 2, • • ■, N. Finally, the APP of aj is given by 

p( a i\ Y) oc p (ai) 7 (aj) 

We now finish the derivations on the BP message passing 
algorithm for SISO MSDD. 


C. Joint Noncoherent Detection and Channel Decoding 

BP message passing algorithm is also widely used as the 
decoding algorithm for many advanced channel codes, such as 
LDPC codes. Turbo codes 0, (15). It is straightforward to 
integrate BP for SISO MSDD with BP for channel decoding 
under the message passing framework, resulting in a BP 
message passing algorithm for joint noncoherent detection and 
channel decoding. In this subsection, we incorporate channel 
encoding/decoding into our framework. The factor graph of the 
overall system that includes the channel encoding constraint 
is shown in Fig. [4] 

The presence of the constraint on a by channel encoding 
introduces loops onto the factor graph. As a consequence, 
the BP message passing algorithm cannot exactly calculate 
these APPs of interest. On the factor graph withe loops, the 
BP message passing algroithm approximates the calculation 
of APPs in an iterative manner — the messages will be 
passed multiple times on some given edges of the factor 
graph Q5). Usually, a good channel code is designed to make 
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Fig. 4. The factor graph for performing the proposed joint MSDD and channel decoding schme, where the factor graph in the dash line box is for performing 
the proposed SISO MSDD, the red arrows denote the intrinsic information sent from SISO MSDD to the channel decoding and the blue arrows denote the 
intrinsic information sent from the channel decoding to SISO MSDD. 


the loops very large. Therefore, in many applications, the 
approximations by iterative BP message passing on factor 
graphs with loops are pretty good ( 13 - 

Given the factor graph of interest, we can design many 
different message-passing schedules. In this work, we adopt 
a serial schedule (3 for the iterative BP message passing 
between the SISO MSDD and the SISO channel decoding. 
The messages exchanged between the MSDD and the channel 
decoding are known as the extrinsic information. In each itera¬ 
tion, given the channel evidences p (Y,; |x^) from observations 
Y and the messages £ (a,.) sent from the channel decoding for 
all i, we performs the BP message passing algorithm for SISO 
MSDD to update messages 7 (a^) as in (17 1 for all i. The 
update of 7 (a^) is still according to (15 1 , (16 1 and (17 1 with 
the only difference that we replace the a priori information 
p(di) in (15 1 and (I 61 with the messages C ( a t)- Then, the 
updated 7 (a,) are treated as the extrinsic information and 
delivered to the channel decoding (denoted by the blue solid 
arrows in Fig. [4|. Then, treating messages 7 (a,;) as the a 
priori information, the channel decoding runs several rounds 
of iterative BP message passing within the subgraph of the 
channel encoding constraint. After that, the channel decoding 
sends back its extrinsic information £ (a*) (denoted by the 
blue dot arrows in Fig. [4]» to the MSDD for the next iteration. 
After several iterations between the MSDD and the channel 
decoding, we terminate the algorithm, and obtain the final 
decoding results about information bits. Finally, we remark 


that the above iterative processing is implemented in digital 
domain as long as we have obtained the correlation samples 
from the AcR receiver. 

D. Alternative SISO MSDD Scheme 

In Section III.B, we establish a hidden Markov chain for 
the signal model induced by the sampling mechanism of 
the proposed AcR; we then develope a BP message passing 
algorithm on the factor graph of the hidden Markov chain 
for fulfilling SISO MSDD. Henceforth, we will refer to this 
hidden Markov chain based MSDD as the M-MSDD scheme. 

In there is another kind of AcR proposed to gen¬ 

erate correlation samples for MSDD. The sampling mechanism 
of the AcR and the following MSDD scheme in G3-G3 is 
in a block-by-block basis. We refer to this block based MSDD 
scheme as the B-MSDD scheme. Originally, the B-MSDD 
scheme proposed in |1()| |.12| only gives hard outputs. It can 
also be modified to become SISO B-MSDD and integrated 
with SISO channel decoding, as investigated in fl8) . In this 
subsection, we interpret the SISO B-MSDD scheme using the 
factor graph and BP message passing algorithm framework. 

The whole packet of N data symbols a*, i = 1, • ■ ■, M is 
divided into U = N/M blocks, where M is the block size. 
The u th block, u = 1, • • •, U, includes M data symbols a*, i = 
(u — 1) M+ 1, • • •, uM. We stack the symbols of the u th block 
into vector x u = [a(„_ 1)M+1 , a (u -i) M + 2 , • ■ ■, a uM \ ■ For B- 
MSDD, the data symbols in the same block will be detected 





























































7 



Fig. 5. The factor graph for performing joint B-MSDD and channel decoding 
with block size M = 3, where i = (u — 1)M + 1, the red arrows denote 
the intrinsic information sent from SISO B-MSDD to the channel decoding 
and the blue arrows denote the intrinsic information sent from the channel 
decoding to SISO B-MSDD. 


jointly. However, different blocks are processed independently. 
For the u th block, the AcR for B-MSDD samples the received 
signal within the observation window with size M +1 symbol 
durations, r (t) for (u — 1) MT S < t < (uM + 1) T s , using 
the following sampling mechanism 



= n 


*=i+1 


N]E g 


°l,J 1 


i = (u — 1) M + 1, • • •, uM 
j = (u - 1) M, ■ ■ ■, i - 1 


08 ) 


where y (t) is the de-spreading signal given in <|5]), n t j is 
the discrete noise component. We stack all the samples Y tJ 
of the u th block into a vector Y„ of length Y^m= t m = 
M (M + l)/2. Then, the observation window is slid down 
M symbol durations for the sampling operations of the next 


M symbols. Based on Y„, B-MSDD proposed in |10)-[12| 
finds hard decisons on the symbols in x u jointly. 


For an SISO B-MSDD, we first construct the probabilistic 
model of the system using factor graphs. The factor graph for 
B-MSDD with M = 3 as an example is shown in Fig. [5] where 
we also incoporate the constraint by the channel encoding. 
Then, we develop a BP message passing algorithm on the 
factor graph. To be integrated with the SISO channel decoding, 
the SISO B-MSDD scheme aims to calculate messages A (a*) 
(denoted by the red arrows in Fig. 0: 


A( Qi) = Y p(y„|x u ) c (aj) (19) 


for * = (u — 1) M + 1, • ■ •, uM, u = 1, • ■ •, U, where I u = 


{(m — 1) M T 1, • • •, uM}, the channel evidence information 



uM i— 1 

oc n n ex p 

i=(u— 1)M+1 j=(u—l)M 

is obtained by using the Gaussian approximation on the 
discrete noise component, ((aj) is the a priori information 
of aj sent from the channel decoding (denoted by the blue 
arrows in Fig. [5]). The messages A (a, ) are treated as the 
extrinsic inforamtion of B-MSDD and sent to the channel 
channel decoding. Then, treating messages A () as the a 
priori information, the channel decoding updates its extrinsic 
information £ ( ai ) and send the updated £ (a^) back to B- 
MSDD. Several iterations between the B-MSDD and the 
channel decoding are performed to ensure the convergence 
of the joint B-MSDD and channel decoding algorithm. 

The SISO B-MSDD scheme and the joint B-MSDD and 
channel decoding scheme are also investigated in fl 8 ) . The 
extrinsic information of B-MSDD in m is computed in 
the logarithm domain. The main contribution of is a 
practical solution to compute the logarithm-domain extrinsic 
information for a large block size M using the list sphere 
decoding | |23) . Here, we just make an interpretation about 
SISO B-MSDD using the framework of the factor graph and 
the BP message passing algorithm. Thus, we compute the 
extrinsic information in the probability domain using the sum- 
product rule as in ( fl9| . The performance of B-MSDD is 
treated as the reference for our M-MSDD. We compare the 
performances of M-MSDD and B-MSDD through simulation 
study in the next section. 

IV. Simulation Results 

In this section, simulations are conducted to validate 
the proposed scheme. In all simulations, the channel are 
generated according to IEEE 802.15.3a CM2 model (24), and 
the channel impulse responses are truncated at T g = 100 ns. 
The used impulse shape uj(t) is the second derivative of a 
Gaussian function. The duration of uj(t') is set as T u = 0.5 
ns. The frame and chip duration are set to Tj = 200 ns 
and T c = 1.0 ns, respectively. Each symbol consists of 
Nf = 10 frames. The TH codes are randomly picked up 
in the interval [0, iV c — 1] where N c = 100. Since now 
T f > T g + (N c — 1 )T C is satisfied, there is no IFI in our 
system. The integration interval of AcR is = T g = 100 ns. 
The bandwidth of the baseband filter employed at the receiver 
is 2 GHz. The signal-to-noise ration (SNR) is defined as 
NfE g /(RN 0 ) 

Test Case 7: 

This case is used to investigate the Gaussian approximation 
on the discrete noise components given in 0 . We get the 
correlation samples from the UWB signal waveforms using 
the sampling mechanism in 0 with M = 3. Then, we 
subtract (iliUi-m+i a *) N f E g f rom Yi,m to get the noise 
component ri l/rn for all m and i. We generate many noise 


















components from 10 4 pakects, each includes N = 1600 data 
symbols (thus, Nf(N + 1) pulses). Finally, we compute the 
empirical noise probability density function (PDF) using 
these simulated noise components. The results are shown in 
Fig. § The theoretical Gaussain PDF with mean zero and 
variance erjj = NfNoE g + WT g N^ /2 are also shown in 
Fig. [6] for comparison. From the results in Fig. [6} we can 
see that the empirical noise PDFs can match the theoretical 
Gaussian PDFs. 

Test Case 2: 

We study whether the used estimation of E g given in ([ 8 ]) 
works well. We compute the mean square error (MSE) of the 
E g estimates. We set that the packet includes N = 1600 data 
symbols, and the block size used to sample the UWB signal 
waveform as in ([2]) is M — 2, 3, 7. For each M and each SNR, 
we average the square errors of E g estimates over 10 4 packets 
to get the MSE. The MSE results are presented in Fig. [7] It 
shows that the values of MSE are small for the median to high 
SNR regime (above 7 dB) and stable for different M. In the 
high SNR regime (above 10 dB), we just see a little difference 
in MSE for different M ; a larger M induces a smaller MSE. 
This is because we can get more samples from lager M and 
we can average out the noises better. Since the MSEs are 
relatively small, we will see later that the use of estimated 
Eg in our scheme just introduce a negligible performance loss. 

Test Case 3: 

We now investigate the bit error rate (BER) performance 
of the proposed M-MSDD scheme for uncoded differential 
UWB-IR systems. Without considering channel codes, we 
employ the M-MSDD to detect the data symbols a* for 
i = 1, 2, • • •, N. After the bidirectional message passing as in 
( fl5| >, ( fl 6 | and ( fi~7j ), the BP message passing algorithm for M- 
MSDD outputs the APP of a, . Then, the hard decision about 
ai is made based upon the APP of a*. Each packet consists of 
N = 1600 data symbols. We evaluate the BER performance 
of the proposed M-MSDD scheme with the perfect E g and 
the estimated E g . As benchmarks, we also evaluate the BER 
performances of the DD (7| scheme and the hard B-MSDD 
1101 |fj2| scheme for differential UWB-IR systems. 

Fig. [ 8 ] presents the BER results. The first point we want to 
study is the impact of E g on the BER performance of our 
M-MSDD. We can observe from Fig. [ 8 ] that the M-MSDD 
schemes with perfect E g and estimated E g nearly have 
the same performance. Thus, we can conclude that the 
estimate of E g by the simple energy estimation method is 
sufficiently effective for the implementation of M-MSDD. 
Then, compared to the DD scheme, our M-MSDD scheme can 
improve the uncoded BER performance by offering detection 
gains (2, 3, 4 dB at the BER of 1(T 6 for M = 2,3, 7) which 
increase with M. This performance trend in improving BER 
by M-MSDD is similar to that by B-MSDD. Since for a fixed 
M, M-MSDD and B-MSDD both have complexity in the 
order of 0(2 M ), we compare the uncoded BER performances 
of M-MSDD and B-MSDD with the same block size M. We 
see that M-MSDD has a better uncoded BER performance 
than B-MSDD — M-MSDD obtains about 1 dB gain at the 


BER of 10 6 for a fixed block size M. 


Test Case 4\ 

We then investigate the BER performances of the joint 
MSDD and channel decoding schemes for coded differential 
UWB-IR systems. The LDPC code |25) with coding rate 
R = 1/2 is employed. Each packet has 800 information bits 
(thus 1600 channel-coded data symbols). In the joint MSDD 
and channel decoding schemes, the used MSDD schemes are 
our SISO M-MSDD proposed in Section III.B and the SISO 
B-MSDD discussed in Section III.D, respectively. We evaluate 
the BER performance of the schemes with the perfect E g 
and the estimated E g . For all simulation results, we perform 
10 iterations between the BP message passing algorithm for 
MSDD and the BP message passing algorithm for LDPC 
channel decoding, and 10 iterations within the BP message 
passing algorithm for LDPC channel decoding. In all simula¬ 
tions, we observe that these numbers of iterations can ensure 
the convergence of the algorithms. 

Fig. [9] presents the coded BER results. The first observation 
is that for the schemes with a large block size M, the estimated 
E g now induces some performance loss in the coded BER. The 
reason is that the coded system with a large M is operating 
at a relatively low SNR regime, where the estimation error of 
E g is large and the decoding performance is sensitive to the 
estimation error. Second, we can see that the proposed joint M- 
MSDD and channel decoding scheme has around 0.4 dB SNR 
gain over the joint B-MSDD and channel decoding scheme for 
a fixed block size M. We believe that these gains are due to the 
more beliefs collected by the BP message passing algorithm 
for M-MSDD. 

In order to show the convergence of the iterative process 
between M-MSDD and the channel decoding, we track the cor¬ 
rect probability of the symbol decisions for a*, i = 1, 2, • • •, N 


through iterations. The results are shown in Fig. 10 The x axis 


represents P c (the channel decoding) , the correct probability 
of the symbol decisions at the output of the channel decoding; 
the y axis represents P c (M — MSDD), the correct probability 
of the symbol decisions at the output of M-MSDD. The 
decisions are made based on the extrinsic information at the 
outputs. The coordinates of each point indicate the correct 
probability of the symbol decisions at the output of M-MSDD 
and the channel decoding at the end of a certain iteration. 
From the results, we find that at a lager SNR operating point, 
the algorithm have a faster convergence. For example, when 
SNR = 8.8 dB, the algorithm with M = 2 has converged 
nearly after 8 iterations between M-MSDD and the channel 
decoding; when SNR = 9.2 dB, the algorithm with M = 2 
has converged nearly after 3 iterations. 


V. Conclusion 

In this paper, we apply BP message passing algorithm to 
propose a joint MSDD and channel decoding scheme for 
noncoherent UWB-IR systems. Specifically, we propose a new 
AcR architecture to transform the received UWB-IR signal 
into discrete samples, whose probabilistic model has a hidden 
Markov chain structure. Using the factor graph representation 
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Fig. 6. Comparison of the empirical noise PDFs to the theoretical Gaussian 
PDFs. The block sampling size is M = 3. 
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Fig. 9. The BER results of the coded system. 



SNR (dB) 

Fig. 7. The MSE results of E g estimates. The used estimation method is 
given in J8j. 
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Fig. 8. The BER results of the uncoded system. 


of this hidden Markov chain model and applying the BP 
message passing algorithm on the factor graph, we derive a 


new SISO MSDD for computing the APPs of the data sym¬ 
bols. The proposed MSDD is a bidirectional message passing 
algorithm, which can make use of all the signal dependences 
throughout the whole packet. Thus, the proposed MSDD has 
a better performance than the block-by-block MSDD scheme 
proposed in previous works. We can also feed the outputs of 
MSDD to the inputs of the BP message passing algorithm for 
channel decoding, and vice verse, in an iterative manner. Our 
simulations indicate that the proposed MSDD scheme has 1 
dB, and 0.4 dB gain over the block-by-block MSDD scheme 
for the uncoded and coded system, respectively. 
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